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CO Abstract. For a detailed comparison of the appearance of cluster of galaxies in X- 

rays and in the optical, we have compiled a comprehensive database of X-ray and optical 

\ properties of a sample of clusters based on the largest available X-ray and optical surveys: 

^ '. the ROSAT AU Sky Survey (RASS) and the Sloan Digital Sky Survey (SDSS). The X-ray 

I ■ galaxy clusters of this RASS-SDSS catalog cover a wide range of masses, from groups 

^ of 10^^-^ Mq to massive clusters of 10^^ Mq in the redshift range from 0.002 to 0.45. 

^ The RASS-SDSS sample comprises all the X-ray selected objects already observed by the 

^ Sloan Digital Sky Survey (114 clusters). For each system we have uniformly determined 
•I— I 

' the X-ray (luminosity in the ROSAT band, bolometric luminosity, center coordinates) and 

$_( 

d ' optical properties (Schechter luminosity function parameters, luminosity, central galaxy 
density, core, total and half-light radii). For a subsample of 53 clusters we have also 
compiled the temperature and the iron abundance from the literature. The total optical 
luminosity can be determined with a typical uncertainty of 20% with a result independent 
of the choice of local or global background subtraction. We searched for parameters which 
provide the best correlation between the X-ray luminosity and the optical properties and 
found that the z band luminosity determined within a cluster aperture of 0.5 Mpc h^Q^ 
provides the best correlation with a scatter of about 60-70%. The scatter decreases to 
less than 40% if the correlation is limited to the bright X-ray clusters. The resulting 
correlation of Lx and Lop in the z and i bands shows a logarithmic slope of 0.38, a 
value not consistent with the assumption of a constant M/L. Consistency is found, 
however, for an increasing M/L with luminosity as suggested by other observations. We 
also investigated the correlation between L^p and the X-ray temperature obtaining the 
same result. 
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1. Introduction 

Cluster of galaxies are the largest well defined building blocks of our Universe. They form 
via gravitational collapse of cosmic matter over a region of several megaparsccs. Cosmic 
baryons, which represent approximately f0-f5% of the mass content of the Universe, 
follow dynamically the dominant dark matter during the collapse. As a resiilt of adiabatic 
compression and of shocks generated by supersonic motions, a thin hot gas permeates 
the cluster gravitational potential. For a typical cluster mass of lO^'^ Mq the intracluster 
gas reaches a temperature of the order of 10^ keV and, thus, radiates optically thin 
thermal bremsstrahlung and line radiation in the X-ray band. In 1978, the launch of 
the first X-ray imaging telescope, the Einstein observatory, began a new era of cluster 
discovery, as clusters proved to be huninous (> 10^2-45 ^^.g^ g-i)^ extended (r ^ 1 — 5 
Mpc) X-ray sources, readily identified in the X-ray sky. Therefore, X-ray observations 
of galaxy clusters provide an efficient and physically motivated method of identification 
of these structures. The X-ray selection is more robust against contamination along the 
line of sight than traditional optical methods since the richest clusters are relatively rare 
and since X-ray emission, which is proportional to the gas density squared, is far more 
sensitive to physical overdcnsities than in the projected number density of galaxies in the 
sky. In fact the existence of diffuse, very hot X-ray emitting gas implies the existence of a 
massive confining dark matter halo. Moreover, selection according to X-ray luminosity is 
also an efficient way to find the highest mass concentrations due to well defined correlation 
between the X-ray luminosity and the total cluster mass (Reiprich & Bhoringher. 2002). 
In addition to allowing the identification of galaxy clusters. X-ray observations provide a 
wealth of information on the intracluster medium itself, e.g. its metal abimdancc, radial 
density distribution and temperature profile. These latter quantities, in turn, can be used 
to reliably estimate the total mass of the system. 

In addition to the hot, diffuse component, baryons are also concentrated in the in- 
dividual galaxies within the cluster. These are best studied through photometric and 
spectroscopic optical surveys, which provide essential information about luminosity, mor- 
phology, stellar population and age. Solid observational evidences indicate a strong inter- 
action between the two baryonic components, as galaxies pollute the intracluster medium 
expelling metals via galactic winds producing the observed metal abundances in clusters 
(De Grandi et al. 2002, Finoguenov et al. 2001). On the other hand, the evolution of galax- 
ies in clusters is influenced by processes due to the hot gas (e.g. gas stripping by ram 
pressure, etc.) as it is by internal processes like star formation, galactic winds, supernovae 
explosions etc., operating inside galaxies themselves (Dressier et al. 1997, Fasano et al. 
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2000, van Dokkum et al. 2000, Lubin et al. 2002, Kelson et al. 1997,2000, Ziegler&Bender 
1997, Gomez et al. 2003). In conclusion, understanding the complex physics at play in 
clusters of galaxies requires combined X-ray and optical observation of a statistically 
significant sample of these objects. 

On the basis of these considerations, we have created a large database of clusters 
of galaxies based on the largest available X-ray and optical surveys: the ROSAT All 
Sky Survey (RASS), the only X-ray all sky survey ever realized using an imaging X- 
ray telescope, and the Sloan Digital Sky Survey (SDSS), which is observing the whole 
Northern Galactic Cap and part of the Southern Galactic Cap in five wide optical bands 
covering the entire optical range. By carefully combining the data of the two surveys 
we have created the RASS-SDSS galaxy cluster catalog. Although two galaxy cluster 
catalogs from the SDSS already exist, the Cut and Enhance Galaxy Cluster Catalog of 
Goto et al. (2002) and the Merged Cluster Catalog of Bahcall et al. (2003, see also Kim 
et al. 2002), we prefered to compile a new cluster catalog by selecting the systems in the 
X-ray band, for which we have reliable X-ray characteristics and for the reasons explained 
above. The X-ray-selected galaxy clusters cover a wide range of masses, from groups of 
1012-5 Mq to massive clusters of 10^^ M© in a redshift range from 0.002 to 0.45. The 
RASS-SDSS sample comprises all the X-ray detected objects already observed in the sky 
region covered by the Sloan Digital Sky Survey. 

One of the first goals is the comparison of the X-ray and the optical appearance of 
the clusters. We want in particular find optical parameters that provide the closest cor- 
relation to the X-ray parameters, such that we can predict within narrow uncertainty 
limits the X-ray luminosity from these optical parameters and vice versa. So far optical 
and X-ray cluster surveys have been conducted independently without much intercom- 
parison. Therefore, the empirical relationship between the X-ray luminosity and optical 
luminosity of clusters is not so well defined, in large part because of the difficulties inher- 
ent in measuring the cluster optical luminosity and in getting a homogeneous set of total 
optical luminosities for a large number of clusters. Edge and Stewart (1991) found that 
the bolomctric X-ray Imninosity of a local sample of X-ray-selectcd clusters correlated 
very roughly with Abell number and somewhat better with the Bahcall galaxy density 
(number of bright galaxies within 0.5 h^^ Mpc; Bahcall 1977,1981) for the small sub- 
sample that had Bahcall galaxy densities. Arnaud et al. (1992) made an heroic effort in 
computing cluster optical luminosities at low redshift from a heterogenous literature. The 
first joint X-ray/optical survey of galaxy clusters was the ROSAT Optical X-ray Survey 
(ROXS, Donahue et al. 2002). They observed 23 ROSAT pointings for a total of 5 square 
degrees in the I band and partially in the V band. The X-ray selection and the optical 
selection of cluster candidates were done independently, with the wavelets algorithm in 
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the former case and with a matched filter algorithm in the latter one. They found X-ray 
and optical coincident detections for 26 galaxy clusters. Donahue et al. (2001) studied 
the relation between the X-ray luminosity and the matched filter parameter Kd, which 
is approximately equivalent to the number of L* galaxies in the system (Postman et al. 
1996). They found a marginally significant correlation between the two quantities with 
a significant scatter. Yee and EUingson (2003) defined a new richness parameter as the 
number of cluster galaxies within some fixed aperture, scaled by a luminosity function 
and a spatial distribution function. They analysed a sample of 15 clusters from CNOCl 
Cluster Redshift Survey, and found a very poor correlation between this parameter and 
other cluster properties such as the X-ray luminosity, temperature and the velocity dis- 
persion. 

In the present paper we describe the properties and the information contained in 
the RASS-SDSS catalog and the resulting correlations between the X-ray and optical 
properties in the sample. In section 2 we explain how the cluster sample has been created 
by X-ray selecting the systems from the available X-ray cluster and group catalogs. In 
section 3 we describe the method for calculating the X-ray cluster properties. We describe 
in section 4 the optical data and in section 5 the data reduction method. We analyse 
and discuss the correlations between the optical luminosity and the X-ray properties in 
section 6. We summarize and discuss the catalog properties and the results in section 7. 

For all derived quantities, we have used Hq = 70 kms~^Mpc~^, Q,m = 0.3 and Oa = 

0.7. 

2. The construction of the sample. 

In order to correlate optical and X-ray galaxy cluster properties it is necessary to have 
a large statistical sample and to cover the whole mass range of the systems considered. 
Since the X-ray observations provide a robust method of identification of galaxy clusters 
and the X-ray luminosity is a good estimator of the system total mass, in principle 
the best approach should bo constructing a cluster catalog of X-ray selected objects in 
a wide range of X-ray luminosity. While for intermediate and high X-ray luminosity 
(mass), several complete catalogs of X-ray selected clusters already exist (NORAS and 
REFLEX), in the low mass range a systematic X-ray survey of groups has still to be 
carried out. As a consequence, it is impossible at the moment constructing a strictly X- 
ray selected cluster sample, which covers a wide range of masses from very poor groups to 
rich clusters. A reasonable compromise in order to fill the low mass range of the spectrum 
is to select all the low mass clusters and groups to date X-ray detected, even if they are 
selected in other wavebands. This compromise is acceptable for our purposes, since we 
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do not want to carry out an X-ray survey of galaxy clusters, but a systematic analysis of 
the correlation between X-ray and optical properties of those systems. 

By following these criteria, the intermediate mass and high mass clusters have been 
selected from three ROSAT based cluster samples: the ROSAT-ESO flux limited X-ray 
cluster sample (REFLEX, Bohringer et al. 2003), the Northern ROSAT All-sky cluster 
sample (NORAS, Bohringer et al. 2000) and the NORAS 2 cluster sample (Retzlaff 
2001) . REFLEX is a complete sample of clusters detected in 13924 deg^ in the southern 
hemisphere down to a flux limit of 3 10~^^ ergs~^cm~^ in the 0.1 — 2.4 keV comprising 
448 clusters. The NORAS galaxy cluster survey contains 495 clusters showing extended 
emission in the RASS in the northern hemisphere with count rates Cx > 0.06 counts 
s~-^ in the 0.1 - 2.4 keV. NORAS 2 is the continuation of the NORAS project and aims 
at a complete survey of X-ray galaxy clusters, in 13598 deg^ of the northern hemisphere, 
down to a flux limit of 2 10~^^ erg s~^cm~^ in the same X-ray band, which gives rise 
to an expected total of about 800 clusters. The samples are based on an MPE internal 
X-ray source catalog extracted with a detection likelihood > 7. 

The low mass clusters and the groups have been selected from two catalogs of X-ray 
detected objects: the ASCA Cluster Catalog (ACC) from Horner et al. (2001) and the 
Group Sample (GS) of Mulchaey et al. 2003. The ACC is a collection of all the clus- 
ters retrieved in the ASCA archive and discovered with different selection strategies. It 
contains measured luminosities, average temperatures, and metal abundances for 273 
clusters and groups. The GS is a heterogeneous collection of 66 ROSAT systems, espe- 
cially optically selected, with velocity dispersions less than 600 km s~^ or an intragroup 
medium temperature less than 2 keV. 

The RASS-SDSS galaxy cluster sample comprises all the X-ray clusters of the cited 
catalogs in the area covered by the Sloan Digital Sky Survey up to February 2003. For 
each X-ray system in the common RASS-SDSS area we found an optical counterpart. The 
sample includes 114 galaxy clusters, 14 of which come from REFLEX, 72 from NORAS 
and NORAS 2, 8 from Mulchaey's groups sample and 20 from ACC. The RASS-SDSS 
galaxy cluster sample, therefore, can not be considered strictly an X-ray 'selected' cluster 
sample, but should be defined more precisely an X-ray 'detected' cluster sample for the 
heterogeneous selection of the low mass range systems. 

3. The X-ray data. 

In order to create a homogeneous catalog of X-ray cluster properties, we have calculated 
all X-ray parameters using only RASS data for all clusters in the sample. The X-ray 
luminosity has been calculated with the growth curve analysis (GCA) method used for 
REFLEX and NORAS 2, based on the RASS3 database (Voges et al. 1999). The GCA 
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method was applied to RASS3 at the position of the clusters. The method allowed for 
a small adjustment of the position to center on the X-ray maximum. Since within the 
GCA aperture the measured flux is underestimated typically by an amount of 7 — 10% 
(Bohringer et al. 2000) the missing flux is estimated by assuming an X-ray luminosity 
scaled cluster model (Bohringer et al. 2001) and corrected for. A first approximate unab- 
sorbed flux is calculated for each X-ray source from the observed count rate, by assuming 
a thermal spectrum with a temperature of 5 keV and a metallicity of 0.3 solar and with- 
out a K-correction. Then an iterative computation uses the redshift and the unabsorbed 
X-ray flux to give a first estimate of the luminosity. With the luminosity-temperature 
relation of Markevitch (1998, without correction for cooling fiows) a better temperature 
estimate is obtained, and the count rate-fiux conversion factor is recomputed including 
the appropriate K-correction for the redshift, resulting in a corrected flux and a new 
X-ray luminosity. The X-ray luminosities as used in this paper are calculated in the 
ROSAT (0.1 — 2.4) keV energy band in the cluster rest-frame for a A cosmology with the 
parameter given above. The GCA also returns for each source many physical parameters 
like improved source position, background brightness, spectral hardness ratio, and KS 
probability for source extent. 

The X-ray bolometric luminosity has been derived from the X-ray luminosity in the 
ROSAT (0.1 — 2.4) keV energy band. A flrst estimation of the cluster temperature is 
calculated by using the Lx — Tx relation of Xue & Wu 2000 in order to estimate the 
appropriate bolometric correction. 

For a subsample of 53 galaxy clusters we have also compiled from the literature the 
ASCA and XMM temperature and iron abundance of the intracluster medium in the 
system. 



4. Optical data 

The optical photometric data were taken from the SDSS (York et al. 2000 and Stoughton 
et al. 2002). The SDSS consists of an imaging survey of tt steradians of the northern sky in 
the flve passbands u, g, r ,i, z, in the entire optical range from the atmospheric ultraviolet 
cutoff in the blue to the sensitivity limit of silicon in the red (Fig. QJ. The survey is 
carried out using a 2.5 m telescope, an imaging mosaic camera with 30 CCDs, two flber- 
fed spectrographs and a 0.5 m telescope for the photometric calibration. The imaging 
survey is taken in drift-scan mode. The imaging data are processed with a photometric 
pipeline (PHOTO) specially written for the SDSS data. 
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4-1- The galaxy sample 

In order to study the optical cluster properties, we have created a complete galaxy sample 
for each cluster, by selecting the galaxies in an area centered on the X-ray source position 
and with radius of 1.5 deg. We used the selection criteria of Yasuda et al. (2001) to 
define our galaxy sample from the photometric catalog produced by PHOTO. We have 
selected only objects flagged with PRIMARY, in order to avoid multiple detections in 
the overlap between adjacent scan lines in two strips of a stripe and between adjacent 
frames. Objects with multiple peaks (parent) are divided by the deblender in diff'erent 
components (children); if the objects can not be deblended, they are additionally flagged 
as NODEBLEND. Only isolated objects, child objects and NODEBLEND flagged objects 
are used in constructing our galaxy sample. The star-galaxy separation is performed in 
each band using an empirical technique, based on the difference between the model and 
the PSF magnitude. An object is classified as galaxy if the model magnitude and the 
Point Spread Function (PSF) magnitude differ by more than 0.145. This method seems 
to be robust for r < 21 mag, which is also the completeness limit of the survey in the 
Northern galactic cap. Since saturated pixels and diffraction spikes can compromise the 
model-fitting algorithm, some stars can be misclassified as galaxies. Therefore we have 
rejected all object with saturated pixels and which are fiagged as BRIGHT. Furthermore, 
we have classified an object as a galaxy only if PHOTO classified it as a galaxy in at least 
two of the three photometric bands g, r, i. After a visual inspection of a sample of galaxies 
with r < 16 mag, Yasuda et al. (2001) concluded that the described selection criteria 
give a sample completeness of 97%, and the same completeness is found for the sample 
at 16 < r < 21 mag after comparison with the Medium Deep Survey catalog (MDS) 
constructed using WFPC2 parallel images from HST. The major reason of missing real 
galaxies is the rejection of galaxies blended with saturated stars, while spurious galaxy 
detections are double stars or shredded galaxies with substructures. 

4-2. SDSS Galaxy photometry 

Since the galaxies do not have sharp edges or a unique surface brightness profile, it 
is nontrivial to define a flux for each object. PHOTO calculates a number of different 
magnitudes for each object: model magnitude, Petrosian magnitudes and PSF magni- 
tudes. The model magnitudes are calculated by fitting de Vaucouleurs and exponential 
model, convolved with the local PSF, to the two dimensional images of the galaxies 
in the r band. The total magnitude are determined from the better fitting of the two 
shape function. Galaxy colors are measured by applying the best fit model of an object 
in the r band to the others bands, measuring the fiux in the same effective aperture. 
However due to a bug of PHOTO, found during the completion of DRl, the model 
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magnitudes are systematically under-estimated by about 0.2 magnitudes for galaxies 
brighter then 20th magnitude, and accordingly the measured radii are systematically too 
large (http://www.sdss.org/DRl/products/catalogs/ind ex.htmll ). This error does not af- 
fect the galaxy colors but makes the model magnitude useless for the determination of 
the galaxy total luminosities. 

The Petrosian flux is defined by 

Fp = 271 I{r)dr, (1) 



where I(r) is the surface brightness profile of the galaxy, and rp is the Petrosian radius 
satisfying the equation: 

„ r-1.25rp I{r)rdr 
r ^ ^^JQ.Srp 7r[(1.25rf.)^-(0.8rf.)^] 

■'1 r) rrp I(r)rdr ' ^ ' 

The Petrosian aperture is set to 2rp, and it encompasses 99% of the galaxy total hght 
in case of an exponential profile and 82% in case of a de Vaucouleurs profile (Blanton et 
al 2001). The Petrosian ratio /i is set to 0.2; at smaller values PHOTO fails to measure 
the Petrosian ratio, since the S/N is too low. For faint objects the effect of the seeing 
on Petrosian magnitude is not negligible. As the galaxy size becomes comparable to the 
seeing disk, the fraction of light measured by the Petrosian quantities approaches the 
fraction for a PSF, about 95%, in which case the flux is reduced for a galaxy with an 
exponential profile and increased for a galaxy with a de Vaucouleurs profile (Strauss 
et al. 2002). Thus the Petrosian magnitudes are the best measure of the total light for 
bright galaxies, but fail to be a good measure for faint objects. In the data analysis of 
this paper we used the Petrosian magnitudes for galaxies brighter then 20 mag and the 
psf magnitudes for objects fainter than 20 mag. 

5. Optical Luminosity from SDSS data 

5.1. Background subtraction 

The total optical luminosity of a cluster has to be calculated after the subtraction of the 
foreground and background galaxy contamination. Since we have used only photometric 
data from the SDSS galaxy catalog, we have no direct information on the cluster galaxy 
memberships. There are two different approaches to overcome this problem. Since galaxy 
clusters show a very well define red sequence in the color magnitude diagram, a galaxy 
color cut could be use to define the cluster membership (Gladders et al. 2000). On 
the other hand the background subtraction can be based on the number counts of the 
projected field galaxies outside the cluster. We chose the latter approach since the former 
method may introduce a bias against bluer cluster members observed to have varying 
fractions due to the Butcher-Oemler effect. 
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We have considered two different approaches to the statistical subtraction of the 
galaxy background. First we have calculated a local background. The M200 — Lx relation 
of Reiprich et al. 2002 was used to compute the r2oo radius (where the cluster mass 
density is 200 times the critical cosmic mass density), as a pragmatic approximation of 
the virial radius. Then, we defined an annulus centered on the cluster X-ray center, with 
an inner radius equal to r2oo + 0.2 deg and a width of 0.5 degree (Fig. In this way 
the galaxy background has been estimated well outside the cluster but still locally. The 
annulus has been then divided in 20 sectors ( analogous to the approach in Bohringcr et 
al. 2001) and those featuring a larger than Scr deviation from the median galaxy density 
are discarded from the further calculation. In this way other clusters close to the target 
or voids are not included in the background correction. We have computed the galaxies 
number counts Nlg{m)dm per bin of magnitude (with a bin width of 0.5 mag) and per 
squared degree in the remaining area of the annulus. The statistical source of error in 
this approach is the Poissonian uncertainty of the counts, given by ^ {Nlg{m)). 

As a second method we have derived a global background correction. The galaxy 
number counts N^g{m)dm was derived from the mean of the magnitude number counts 
determined in five different SDSS sky regions, each with an area of 30 deg^ (FigOl- 
The source of uncertainty in this second case is systematic and originates the presence 
of large-scale clustering within the galaxy sample, while the Poissonian error of the 
galaxy counts is small due to the large area involved. We have estimated this error as 
the standard deviation of the mean global number counts, a^g{m), in the comparison 
of the five areas. In order to take into account this systematic source of error also for 
the the local background , we have estimated the background number counts error as 
o'bg(™) = TOaa;(y^(A^^g(m)), cr^g(»7i)) (Lumdsen et al. 1997) for all the derived quantities. 

After the background subtraction we found that the signal to noise in the u band was 
too low to be useful, and performed our analysis on the 4 remaining Sloan photometric 
bands g,r,i,z. 

5.2. Luminosity Function 

For each cluster and in all photometric bands we have assumed that the distribution 
function of galaxies in magnitude can be described by the Schechter Luminosity function 
(LF): 

(t){m)dm = 0.4 ln[lQ) N^iu <j)* 10-°-* (m-m*) (a-l) ^xp{-lQ-° '^ (m-m*)) (3) 

In the equation Ndu is the number of cluster galaxies and was computed as the difference 
between the total number of galaxies in the cluster region and the expected number 
of interlopers, estimated from the local (global) background galaxy density. 0* is the 
normalization of the Schechter Luminosity function, given by the inverse of the integral 
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of the LF over the considered magnitude range. To determine the remaining parameters 
M* and a we fitted the Schechter LF to the data with a Maximum Likehhood Method 
(MLM, Sarazin f980). Since we have no information about the cluster membership, we 
have considered the observed galaxy magnitude distribution in the cluster region as the 
sum of the Schechter Luminosity function {(j){m)) plus a background contribution (6(m)): 
0(m)+6(m) 

= \r ■ (4) 

$(m) is normalized to unity when integrated over the considered range of magnitudes. 

In order to perform a ML analysis, the background contribution has to be specified at 
any magnitude. To estimate the h{m), one can try to fit the background number counts 
by a model. While the behavior of the N{m) — m relation is well know at the bright 
end (in the logNijn) — logm is a line with slope 0.6, Figl^l , it is not well understood 
at the faint end (Yasuda et al. 2001, Lumdsen et al. 1997). Therefore instead of using 
a specific functional form, we simply have used a spline to interpolate the background 
galaxy number counts and estimated b{m) at any magnitude. 

The probability that the assumed distribution gives a galaxy at the magnitudes ruk 
is thus $(mfc). Therefore, if the observed galaxies are statistically independent, the com- 
bined probability that the assumed distribution gives the observed galaxies at the mag- 
nitude mfc (with k = l,7i) is: 

n 

L=[]$(mfc) (5) 

K=l 

The best-fit parameters are those that maximize the likelihood L. In practice we have 
minimized the log-likelihood -21n(L). This minimization was performed with the CERN's 
software package MINUIT. We used the variable metric method MIGRAD (Flechter 
1970) for the convergence at the minimum and the MINOS routine to estimate the error 
parameters in case of non linearities. We also have placed constraints on the values of 
m* and a that the fitting routine can accept, to avoid being trapped in a false minimum 
(M* in the range between -18 and -26 mag and a between and -2.5, Lumdsen et 
al. 1997). Fig. 0] gives an example of the LF derived with MLM in each of the Sloan 
photometric bands for a cluster. Figures and El show the comparison between the fit 
parameters calculated with different backgrounds h{m) (local and global). Fig. [3 shows 
the correlation between the fit parameters M* and a. 

The great advantage of the MLM is that the method does not require the data to be 
binned and does not depend on the bin size, but uses all the available information . On the 
other hand, the MLM provides no information about the goodness of the fit. Therefore, 
we have performed a statistical test. Since the routine procedure uses a unbinned set of 
data to perform the fit, in principle the Kolmogorov Smirnov test should be applicable to 
our case. Nevertheless since the KS probability is not easy to interpret, we have applied 
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a test, by comparing the background subtracted magnitude number counts of the 
cluster with the Schechter luminosity function fitted to the data. Figure |H1 shows the 
distribution of the reduced hi the cluster sample. Almost 90% of the fitted LFs are a 
good fit to the data having a reduced < 1. 

5.3. The total luminosity 

In order to calculate the total cluster luminosity, we have calculated first the absolute 
magnitude 



where is the luminosity distance, A is the Galactic extinction and K{z) is the K- 
correction. We deredden the Petrosian and model magnitudes of galaxies using the Galac- 
tic map of Schlegel et al. (1998) in each photometric band. We used the K-correction 
suppHed by Fukugita, Shimasaku, & Ichikawa (1995) for elliptical galaxies, assuming that 
the main population of our clusters are the old elliptical galaxies at the cluster redshift. 
The transformation from absolute magnitudes to absolute luminosity in units of solar 
luminosities is performed by using the solar absolute magnitude obtained from the color 
transformation equation from the Johnson-Morgan-Cousins system to the SDSS system 
of Fukugita et al. (1996). We have calculated the optical luminosity of each cluster with 
two different methods. First, we have estimated L by using the (background corrected) 
magnitude number counts of the cluster galaxies with the following prescription: 



The sum on right side is performed over all the N magnitude bins with galaxy number 
Ni{m) and mean luminosity li{m). The integral is an incompleteness correction due to the 
completeness limit of the galaxy sample at = 21 mag in the five Sloan photometric 
bands. (/)(m) is the individual Schechter luminosity function fitted to the galaxy sample of 
each cluster. The incompleteness correction is of the order of 5-10% in the whole cluster 
sample, as showed in the Fig.|^ This means that the galaxies below the magnitude limit 
do not give a significant contribution to the total optical luminosity. Therefore the most 
important source of error is due to the contribution of the background galaxy number 
counts. The uncertainty in each bin of magnitude is given by the Poissonian error of 
the bin counts ( yiV^^^Jro) , with i — 1,...,N) and the background subtraction in each 
magnitude bin {o'lg^m), with i = 1, ...,7V, see previous section for the definition). Since 
the galaxy counts in the bins are independent, the error in the luminosity is given by: 



M = m - 25 - 5logio{DL/lMpc) -A- K{z) 



(6) 




(7) 



N 




(8) 
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Fig. ^1 shows the comparison between the luminosity calculated from the local back- 
ground corrected and global background corrected magnitude number counts. The dif- 
ference between both methods is much smaller than the statistical error. 

In the second case we have taken advantage of the individual fitted luminosity func- 
tion. If MIGRAD has converged successfully to the minimum, we calculate the total 
luminosity as 

L*N,i^(b*r{2 + a) (9) 

where L* and a are the fit parameters estimated from the data and N^u is the total 
number of cluster galaxies. There are different sources of errors in this calculation. The 
major source of error comes from the background which affects both the number of cluster 
galaxies N^u and the result of the fitting procedure. A second kind of error is due to the 
uncertainty of the fit parameters. Since all these errors are not independent, we can not 
treat their contributions separately. Therefore the luminosity errors were calculated by 
varying the fit parameter values, M* and a, along their 68% confidence level error ellipse 
and using the upper and lower bound of the quoted background number counts (b(rn)) 
ranges. The statistical luminosity error range was then defined between the minimum 
and maximum luminosity. With this method we can take into account statistical and 
systematic errors due to the background and their effects on the fit parameters as well. 
Note that a simple error propagation applied to the equation would underestimate 
the error in the luminosity, since it would not take into account the error of the galaxy 
background. 

Fig. ^3 shows the comparison of the two optical luminosities (fit-based and count- 
based), which are consistent within the errors. For 70% of the clusters in the sample the 
count-based luminosity is systematically brighter than the the fit-based one, as showed in 
the Fig.^] In the former case, indeed, the method includes in the calculation of Lop the 
Bright Cluster Galaxies (BCG) which are usually excluded by the Schechter luminosity 
function. The error bars in the fit-based luminosity are larger than in the count-based 
ones. In fact in the former case there are two main sources of error: the uncertainty due 
to the galaxy background subtraction and the statistical errors in the fit parameters of 
the luminosity function. In the latter case. Instead, only the subtraction of the galactic 
background plays a crucial role. The mean error in the fit-based luminosity is around 
30%, while it is around 20% with the count-based method. 

The great advantages of the count-based optical luminosity is that it can be easily 
computed, if the cluster S/N is high enough. On the other hand, the fit-based luminosity 
depends on the success of the fitting procedure. Therefore, it is sensitive not only to 
the S/N but also to the chosen model and the goodness of the fit. The uncertainty in 
the count-based method is smaller than in the fit-based method. Moreover, while the 
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count-based method provides the optical luminosity for any system and at any cluster 
aperture, the number of failures in the fitting procedure is an increasing function of the 
cluster aperture. In fact the fit-based method fails to fit the data for 15% of the clusters 
at 0.35 Mpc hyQ^ up to 35% at 2.0 Mpc hyQ^. In consequence, the count-based Lop has 
to be preferred to the fit-based one in the study of the correlation between optical and 
X-ray properties. The count-based method also reflects what we actually observe. 

On the basis of this analysis we can conclude that the behavior of the optical lumi- 
nosities calculated with different background subtraction is stable for variant approaches 
and the main source of errors is due to the necessary background subtraction fFig. llOji . 
Moreover, since the two different methods (count-based and fit-based) give consistent 
results, our measure of Lop seems to be a good estimation of the cluster total optical 
luminosity. 



5. 4- The optical structure parameters 

In order to study the spatial distribution of galaxies in cluster, we have analysed the 
projected radial galaxy distributions of each cluster in the sample. The analysis is per- 
formed in the g,r i and z bands. As for the luminosity functions, we have used a Maximum 
Likelihood method to fit a King profile to the data 

^W= (l+[1)2)^ +"''- (10) 

In eauation llOl erg is the central galaxy density, r^ the core radius, (3 the profile exponent, 
and cTfc the background density. P{r) has to be normalized through: 

P{r)d{^r^) = Ntot (11) 

where A is the relevant cluster area and Ntot is the total number of galaxies within that 
area. In agreement with the section 5.2 the the Likelihood is given by: 

n 

i = n ^(^-) (12) 

where r^ is the projected galaxy distance from the X-ray center. We regarded (3, rc and 
(Tb as fitting parameters, while ctq is a dependent variable and its value is derived from the 
likelihood normalization, eq. (|lll) . The fitting method worked successfully in average for 
95% of the sample in any photometric band; it failed for groups, where the overdensity 
in comparison to the background density is too low to fit a profile. As shown in the fig. 
1131 there are no correlations between the parameters, ctq, rc and /3, with the background 
density at- Furthermore the histogram of the f3 values in the same figure shows that 
the mean value of the profile exponent is around 0.8 with a very large dispersion of 0.5 
around the peak. 
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We have estimated from the King profile the physical size of each cluster, rtot- We have 
assumed that this quantity is given by the radial distance from the X-ray center, where 
the galaxy number density of the cluster becomes n times the error of the background 
galaxy density ( the cross in the Fig. I14() . To search for the best value of n, we have 
estimated the total radius with different values of ?i (n = 1,...,5) and calculated the 
total optical luminosity within that radius; n was then fixed to 3, since the differences 
in the luminosities calculated within different total radii are smaller than the luminosity 
uncertainties due to background subtraction. 

We assumed that the cluster total optical luminosity in a band is the luminosity 
calculated within rtot estimated in the given filter. To calculate then the half-light radius 
, which encircles half of the total cluster luminosity, we have estimated in each band 
the luminosity of the cluster within 20 radii from the X-ray center to the total radius. 
We have then interpolated in each filter the 20 luminosities, to find the radius which 
corresponds to half of the total cluster luminosity. 

5.5. The catalog 

In the following we present the catalog of the 114 RASS-SDSS galaxy clusters. Tables 
1-3 list all the X-ray and optical properties of the sample computed as explained in the 
previous sections. The examples tables show the results for the first 35 clusters in the 
sample. The complete tables are given in electronic form. 

Table give the X-ray properties of the cluster derived from the ROSAT data. 
Col. (1) and (2) contain the ROSAT and the alternative cluster name respectively. Col. 
(3) and (4) contain the equatorial coordinates of the X-ray cluster center used for the 
regional selection for the epoch J2000 in decimal degrees. Col. (5) contains the heliocentric 
cluster redshift. Col. (6) presents the fiux in the energy band range 0.1-2.4 keV in units 
of 10~^^ ergs s~^ cm~^. Col. (7) gives the corrected fiux for a temperature derived from 
Lx — T relation including the K-correction for an assumed cluster temperature of 5 kcV. 
Col. (8) contains the relative 1 cr Poissonian error of the count rate, the fiux and the 
luminosity in percent. Col (9) gives the luminosity in the energy range 0.1-2.4 keV in 
units of hyQ^10'*'*erg s~^. Col. (10) contains the count rate in units of counts s~^. Col. 
(11) gives the outer radius within which the fiux and the luminosity are estimated, in 
units of arc minutes. 

Tables (jJl provides the optical parameters of the luminosity function and the lu- 
minosities of each cluster calculated by using the local galaxy background. Listed are 
results in the r band. All the quantities are calculated within a fixed cluster aperture of 
1.0 Mpc hfo^ Col. (1) presents the ROSAT name of the cluster. Col. (2) and (3) show 
the resulting fit parameters of a Schechtcr luminosity hmction. a is the slope of the LF, 
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while M* is the magnitude knee of the distribution. CoL (4) gives the galaxy number 
density within the cluster region selected to perform the fit, in units of deg^^. Col. (5) 
shows the reduced of the fitted luminosity function. Col (6) provides the cluster optical 
luminosity in units of 10^^ Lq, calculated on the basis of the fitted luminosity function. 
Col (7) lists the cluster optical luminosity in units of 10^^ i©, calculated on the base of 
the cluster magnitude number counts. All the errors in the table are at 68% confidence 
level. The catalog contains the extended version of this table. Similar tables exist, listing 
the parameters and the luminosities relative to 22 different cluster apertures: 20 fixed 
apertures ranging from 0.05 to 4.0 Mpc h^Q^, 2 variable apertures as the core radius and 
the half-light radius. All the data are provided in each of the 4 Sloan photometric bands 
g,r,i and z, and for both local and global galaxy background correction. 

Table Q provides the information concerning to the radial distribution of the pro- 
jected galaxy density in the region of the cluster. Col. (1) presents the ROSAT name of 
the cluster. Col. (2) gives the cluster central galaxy number density in units of deg~^ (tro 
is not a fit parameter, therefore the error is not provided). Col. (3) lists the background 
galaxy number density around the cluster in units of deg~^. Col. (4) shows the core 
radius estimated from the fit, in units of Mpc. Col. (5) provides the cluster total radius, 
extrapolated from the King profile, in units of Mpc. Col. (6) gives the half-light radius 
in units of Mpc. All the errors in the table are at 68% confidence level. 

The full set of extended tables are available in electronic form. 

6. Correlating X-ray and optical properties 

For a cluster in which mass traces optical light (M/Lopt is constant), the gas is in hy- 
rostatic equilibrium (T cx M^/^), and Lx oc (Xue & Wu 2000), we expect the X-ray 
bolometric luminosity to be related to the optical luminosity as Lop oc i*^^ and to the 
intracluster medium temperature as Lop (x T^^. 

We have now an optimal data base to test these scaling relations. In a first step we 
look for those optical parameters which are best suited for a correlation analysis and 
apply these in a second step to the test of the scaling relations. 

In this section we show that tight correlations exist between the total optical cluster 
luminosity and the X-ray cluster properties as the X-ray luminosity and the intracluster 
medium temperature. 

To search for the best correlation between optical and X-ray properties and to opti- 
mally predict for example the X-ray luminosity from the optical appearance, we are in- 
terested in an optical characteristic, which shows a minimum scatter in the X-ray/optical 
correlation. Therefore, we performe a correlation using 4 of the 5 SDSS optical band, g, 
r, i and z, to find out which filter should be used in the prediction. The u band was not 
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Table 1. Example of the table containing the X-ray cluster properties of the whole 
cluster sample. The table show the results for the first 35 cluster in the sample. 
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0, 


.69 


2, 


.417 


± 


0, 


.076 


2, 


.434 


± 


0, 


.717 


RXCJ1022.5+5006 


-0.65 ± 0.72 


-20.52 + 0.65 


5544 


0, 


.21 


2, 


.265 


± 


0, 


.032 


2, 


.348 


± 


0, 


.981 


RXCJ1023.6+0411 


-0.75 ± 0.40 


-20.58 ± 0.73 


2442 


0. 


.60 
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.079 


± 


0, 


.036 


1, 


.309 


+ 


0. 


.492 


RXCJ1023.6+4908 


-1.31 + 0.13 


-21.80 + 0.65 


1227 


1. 


.06 


1. 
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± 


0, 


.010 


1, 
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+ 


0. 


.358 


RXCJ1053.7+5452 


-1.52 + 0.21 
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1. 
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± 


0, 


.203 


2, 
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+ 


0. 


.750 


RXCJ1058.4+5647 


-1.49 + 0.09 


-21.98 + 0.64 


2557 


0. 


.63 


1. 


.338 


± 


0, 


.011 


1, 


.338 


+ 


0. 


.409 



Table 2. Example of the table containing optical properties of the whole cluster sample 
in the r band. The table show the results for the first 35 clusters in the sample. 
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73049 


4563 ± 242 


0, 


.006 


± 


0, 


.003 


2, 


.741 


1, 


.516 


± 


0. 


.178 


RXCJ1023.6+0411 


14797 


4129 ± 62 


0. 


.049 


lb 


0, 


.006 


1, 


.907 


0. 


.337 


± 


0. 


.036 


RXCJ1023.6+4908 


3003 


5215 ± 53 


0. 


.175 


lb 
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.027 


0, 


.712 


0. 
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0. 


.066 


RXCJ1053.7+5452 


22426 


3428 ± 265 
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lb 


0, 
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2, 
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± 


0. 
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5110 ± 135 
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lb 
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.269 
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± 
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.088 



Table 3. Example of the table containing the optical properties of the whole cluster 
sample in the r band. The table shows the results for the first 35 cluster of the sample. 
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used since the cluster signal to noise in that band is too low to calculate the cluster total 
luminosity. We used a fixed aperture to calculate the optical luminosities for all the clus- 
ters, to make no a priori assumption about the cluster size. Moreover, to check whether 
the scatter in the correlation depends on the cluster aperture, we did the same analysis 
several time by using optical luminosities calculated within different radii, ranging from 
0.05 to 4 Mpc h^Q from the X-ray center. To quantify the Lop — Lx and the Lop — 
relations, a linear regression in log-log space was performed by using two methods for the 
fitting: a numerical orthogonal distance regression method (ODRPACK) and the bisector 
method (Akritas & Bershady 1996). The fits are performed using the form 



where Px is the X-ray property, and the errors of each variable are transformed in to 
log space as /S.log{x) — log{e){x^ — x^)/{2x), where and denote the upper and 
lower boundary of the quantity's error range, respectively. To exclude the outliers in the 
fitting procedure we apply a a clipping method. After a first fit all the points featuring a 
larger than 3(7 deviation from the relation, were excluded and the fitting procedure was 
repeated (see discussion below) . 

Figs. 1 151 and 1 161 show the scatter of the Lop — Lx and the Lop — relations, respec- 
tively, versus the cluster aperture, used to calculate the optical luminosity. The scatter 
in the plot is the orthogonal scatter estimated from the best fit given by ODRPACK. 
In any photometric band the scatter has a clear dependence on the cluster aperture by 
showing a region of minimum on the very center of the cluster, between 0.2 Mpc h^Q^ 
and 0.8 Mpc hyQ^, and by increasing at larger radii. The source of the scatter is partially 
due to the method used to calculate the optical luminosity. Our method is simply based 
on the overdensity of the cluster with respect of the galaxy background. In fact if the 
central region in which Lop is measured is to small ( 0.05 Mpc h^Q^ in Figs. 1151 and ll6|l . 
the value of the galaxy density is low and the measurement becomes more uncertain. 
On the other hand, at larger radii the density contrast between cluster and background 
decreases progressively. Instead, within a cluster aperture between 0.2 Mpc Kj^ and 0.8 
Mpc h^Q^ , the optical luminosity of both groups and massive clusters can be easily mea- 
sured. In fact in both cases the radial aperture is small enough to show an high density 
contrast and therefore an high cluster signal to noise, and still large to enclose enough 
galaxies for the luminosity calculation. 

After a more accurate analysis, we noted also that the low luminosity systems (both 
in the optical and in the X-ray band) are the main source of scatter at any cluster 
aperture. This could be due to different reasons. From the technical point of view the 
groups have a lower surface density contrast, and this causes problems in calculating the 



log{Lop/LQ) 



alog{Px)+f} 



(13) 
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optical luminosity with a method based on the overdensity contrast. Moreover, the low 
mass systems could have a larger scatter in the optical and X-ray properties. 

Furthermore, the galaxy groups could be responsible for the behavior of the scatter 
shown in fig. 1151 In fact at large cluster apertures the galaxy density contrast can be 
very low for the small systems and still very high for the massive and larger clusters. 
The large error introduced by the low density contrast in the calculation of the optical 
luminosity of galaxy groups could be the source of scatter able to explain the increment of 
the scatter at larger apertures. In order to study in more details the nature of the scatter 
of our correlation and in order to investigate the role of the less luminous systems, we 
carried out the analysis explained above with the low mass systems removed. We limited 
the analysis to the subsample of the X-ray selected REFLEX-NORAS clusters, which 
occupy the intermediate and high luminosity region. Fig. 1171 shows the behavior of the 
scatter as a function of the cluster aperture in this second analysis. After removing the 
low mass systems the scatter decreases by about 30% for any aperture (from 0.3 dex 
to 0.2 dex in the region of minimum scatter and from 0.6 dex in average to 0.4 dex 
at larger apertures). Nevertheless the behavior of the scatter at increasing aperture is 
absolutely the same observed in the analysis carried out on the overall RASS-SDSS galaxy 
cluster sample. This means that groups are responsible for part of the scatter but can 
not explain the existence of the region of minimum scatter between 0.2 and 0.8 Mpc h^Q ■ 
A possible explanation for the behavior of the scatter at different cluster apertures could 
be the cluster compactness. As Lx depends not only on the cluster mass but also on the 
compactness of the cluster, also the optical luminosity should reflect somehow the cluster 
properties, mass and concentration. Thus, there should be an optimal aperture radius 
within which to measure Lop. We found that in all photometric bands, the minimum 
scatter is around 0.5 Mpc h^Q- 

Figg . ni^ and 1191 show the Lgp — Lx and Lop — relation respectively, at the radius 
of minimum scatter, 0.5 Mpc h^Q^, in the z band. In fact the i and the z bands have a 
slightly smaller scatter than in the other optical bands at any radius. Both plots show, 
as an outliers, the cluster RXCJ0845. 3-1-4430 , which features a deviation larger than 
3(7 from the best fit. The system is a nearby group with a density contrast too low to 
estimate the optical luminosity reliably. The X-ray luminosity of this system has a 100% 
error respectively. In the Lop — T^ there is another outlier: the cluster RXCJ1629. 6-1-4049. 
This system is not a source of scatter in the Lop — Lx relation and the error in the X-ray 
and optical luminosities is 8% and 35% respectively. This can suggests that the optical 
luminosity is well measured, and it questions the estimate of the temperature. In fact, 
the X-ray luminosity of RXCJ1629. 6-1-4049 is 2.78x10^^ ergss"\ and the temperature. 
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estimated from Horner et al. (2001), is 1 keV, while the Lx — Tx relation of Xue & Wu, 
(2000), predicts at least a Tx of 4.3 keV at that Lx- 

With the a clipping method those clusters were excluded from the estimation of the 
best fit. The best fit parameters of the orthogonal and bisector methods, in the i and 
z photometric bands, arc shown in the Table 01 and for the Lop — Lx relation for 
the ROSAT X-ray luminosity and the bolometric X-ray luminosity, respectively. Table 
El shows the same results for the Lop — relation. Table [3 provides the results for the 
Lop — Lx relation in i and z band for the subsample of X-ray selected REFLEX-NORAS 
clusters. The tables show also the estimated orthogonal scatter and the estimated scatter 
in both the variables. The best fit in the z band for the Lop — Lx and the Lop — Tx 
relations at the radius of minimum scatter for all the RASS-SDSS galaxy cluster sample 
are respectively: 

Lop/Lq = 10"™±"02L;f (i?O5Ar)"-45±0"3 (14) 

Lop/Lq = 10"-^5±0.02^^(^^;)0.38±0.02 (^5) 
Lop/Lq = I0ll-42±0.06j,^.12±0.08 (^g) 

The value of the exponent in the power law for the Lopt — Lx{Bol) relation, is around 
0.38 in the region of minimum, as indicated in tabled The values are not consistent 
within the errors with the value of 0.5 predicted under the assumption of hydrostatic 
equilibrium and constant mass to light ratio. The same conclusion can be derived for 
the Lopt — Tx relation and from the Lop — Lx for the subsample of X-ray selected 
REFLEX-NORAS clusters. A simple reason of the disagreement could be due to the 
assumption of a constant mass to light ratio. In fact, Girardi et al. (2002) analysed in 
detail the mass to light ratio in the B band of a sample of 294 clusters and groups , 
finding M/L cx l0.33±o.03^ rpj^^ ^^^^ results was found by Lin et al. (2003) in the K 
band. Thus if we consider this dependence of M/L from the optical luminosity with the 
assumptions of hydrostatic equilibrium, the new expected relation between the optical 
luminosity and the X-ray luminosity and temperature are Lop oc L^^ and Lop oc T}^^'*^ 
respectively, which are in good agreement with our results. 

The scatter in the Lop ~ L^ relation for the aperture with the best correlation ( 
0.5 Mpc h^Q^), in the Lop variable is 0.20, and the scatter in the Lx variable is 0.22 
in the correlations obtained in i and the z bands as shown in the Table 01 Therefore, 
by calculating the total cluster luminosity in the central part of the system, one can 
use the z or z band to predict the X-ray luminosity from the optical data with a mean 
error of 60%. In the same way the optical hmiinosity can be derived from Lx with the 
same uncertainty. As indicated in the table the uncertainty in the prediction of the two 
variables decreases to less than 40% if the correlation is limited to the REFLEX-NORAS 
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Lop — Lx (ROSAT) relation in the I band 



Orthogonal method 



Bisector method 



r 


a 


P 


a 






a 


/3 


a 


CTLx 




0.2 


0.50 ± 0.04 


11.41 ±0.03 


0.43 


0.32 


0.30 


0.32 ±0.01 


11.44 ±0.03 


0.36 


0.23 


0.28 


0.3 


0.46 ±0.03 


11.70 ±0.02 


0.31 


0.23 


0.21 


0.36 ±0.01 


11.71 ± 0.02 


0.29 


0.18 


0.23 


0.5 


0.47 ±0.03 


11.83 ±0.02 


0.30 


0.22 


0.20 


0.38 ±0.01 


11.85 ± 0.02 


0.27 


0.18 


0.20 


0.7 


0.56 ±0.03 


11.94 ±0.03 


0.36 


0.26 


0.25 


0.47 ±0.01 


11.96 ± 0.02 


0.32 


0.22 


0.24 


0.8 


0.61 ±0.04 


11.99 ±0.03 


0.40 


0.29 


0.28 


0.51 ±0.01 


12.01 ± 0.03 


0.36 


0.24 


0.26 


1.0 


0.65 ± 0.04 


12.05 ±0.03 


0.44 


0.31 


0.31 


0.54 ±0.01 


12.07 ±0.03 


0.39 


0.27 


0.28 






Lop 


~Lx 


(ROSAT) relation in the Z band 








Orthogonal method 


Bisector method 


r 


Q 


p 


a 


aLx 




a 




a 


aLx 


^Lap 


0.2 


0.50 ± 0.04 


11.50 ±0.03 


0.43 


0.31 


0.31 


0.31 ±0.01 


11.54 ±0.03 


0.35 


0.23 


0.26 


0.3 


0.45 ± 0.03 


11.79 ±0.02 


0.31 


0.23 


0.21 


0.37 ±0.01 


11.80 ± 0.02 


0.29 


0.18 


0.23 


0.5 


0.47 ±0.03 


11.92 ±0.02 


0.30 


0.22 


0.20 


0.39 ±0.01 


11.93 ± 0.02 


0.28 


0.18 


0.21 


0.7 


0.56 ±0.03 


12.03 ±0.03 


0.36 


0.26 


0.25 


0.47 ±0.01 


12.04 ±0.02 


0.33 


0.22 


0.24 


0.8 


0.60 ± 0.04 


12.08 ±0.03 


0.40 


0.29 


0.27 


0.51 ±0.01 


12.10 ± 0.03 


0.36 


0.24 


0.26 


1.0 


0.65 ± 0.04 


12.14 ±0.03 


0.44 


0.31 


0.31 


0.54 ±0.01 


12.15 ± 0.03 


0.39 


0.27 


0.28 



Table 4. The table presents the results of the best fit for the i and z band in the region 
of the minimum scatter, 0.2 < r < 1.0 Mpc h^Q^ for the Lop — Lx (ROSAT) relation. 
We show the results for the two methods applied: the orthogonal distance regression 
(ODRPACK) and the bisector method, which is the line bisecting the best fit results of 
the vertical and horizontal minimization. The a and (3 parameters are given for several 
apertures with the 95% confidence errors. The orthogonal scatter and the scatters in 
log(Ljs: (0.1-2.4 keV) and in log(Lop) to the best fit line are given by a, (Tl^ and (JLap 
respectively. 



cluster subsample. Table shows that analogous results are obtained for the Lop — 
relation. 

Since the observational uncertainties in the optical and in the X-ray luminosity are 
about 20%, the scatters of 60% of the overall sample in both relations and of 40% in the 
the REFLEX-NORAS cluster subsample for the Lop — L^ relation should be intrinsic. 

7. Summary and Conclusions 

We created a database of clusters of galaxies based on the largest available X-ray and 
optical surveys: the ROSAT AU Sky Survey (RASS) and the Sloan Digital Sky Survey 
(SDSS). The RASS-SDSS galaxy cluster catalog is the first catalog which combines X- 
ray and optical data for a large number (114) of galaxy clusters. The systems are X-ray 
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Lx {Bolometric) relation in the I band 



Orthogonal method 



Bisector method 
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/3 
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/3 


a 






0.2 


0.38 ±0.03 


11.25 ±0.04 


0.41 


0.31 


0.27 


0.27 ±0.01 


11.31 ±0.03 


0.37 


0.24 


0.28 


0.3 


0.36 ±0.02 


11.54 ±0.03 


0.31 


0.23 


0.20 


0.30 ±0.01 


11.57 ±0.02 


0.29 


0.18 


0.23 


0.5 


0.38 ±0.02 


11.66 ±0.03 


0.30 


0.22 


0.20 


0.32 ±0.01 


11.70 ± 0.02 


0.27 
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Table 5. The table presents the results of the best fit for the i and z band in the region 
of the minimum scatter, 0.2 <r< 1.0 Mpc h^Q^ for the Lgp — Lx {Bolometric) relation. 
We show the results for the two methods applied: the orthogonal distance regression 
(ODRPACK) and the bisector method, which is the line bisecting the best fit results of 
the vertical and horizontal minimization. The a and f3 parameters are given for several 
apertures with the 95% confidence errors. The orthogonal scatter and the scatters in 
log(Tx-) and in log(Lop) to the best fit line are given by a, axx and aLop respectively. 



selected, ranging from groups of 10^^ ''' Mq to massive clusters of 10^'"' Mq in a rcdshift 
range from 0.002 to 0.45. The X-ray (luminosity in the ROSAT band, bolometric lumi- 
nosity, redshift, center coordinates) and optical properties (Schcchtcr luminosity fmiction 
parameters, luminosity, central galaxy density, core, total and half-light radii) are com- 
puted in a uniform and accurate way. The catalog contains also temperature and iron 
abundance for a subsample of 53 clusters from the Asca Cluster Catalog and the Group 
Sample. The resulting RASS-SDSS galaxy cluster catalog, then, constitutes an important 
database to study the properties of galaxy clusters and in particular the relation of the 
galaxy population seen in the optical to the properties of the X-ray luminous ICM. 

The first investigations reported have shown a tight correlation between the X-ray and 
optical properties, when the choice of the measurement aperture for the optical luminosity 
and the optical wavelength band are optimized. We found that the optical luminosity 
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Table 6. The table presents the results of the best fit for the i and z band in the region 
of the minimum scatter, 0.2 < r < 1.0 Mpc h^Q^. We show the results for the two methods 
applied: the orthogonal distance regression (ODRPACK) and the bisector method, which 
is the line bisecting the best fit results of the vertical and horizontal minimization. The 
a and /3 parameters are given for several apertures with the 95% confidence errors. The 
orthogonal scatter and the scatters in log(Tx-) and in \og{Lop) to the best fit line are 
given by ct, cttx and aL„f, respectively. 



calculated in the i and in the z band correlates better with the X-ray luminosity and 
the ICM temperature, in comparison to the other Sloan photometric bands. Thus the 
red optical bands, which are more sensitive to the light of the old stellar population and 
therefore to the stellar mass of cluster galaxies, have tight correlations with the X-ray 
properties of the systems. 

Moreover, wo foimd that the scatter in the Lop — Lx and Lgp ~ Tx relations can be 
minimized if the optical luminosity is measured within a cluster aperture between 0.2-0.8 
Mpc h^Q^, with an absolute minimum of the scatter at 0.5 Mpc hjQ . The best aperture 
in the central part of the cluster for the measurement of the optical luminosity is due to 
the fact that it is a good compromise to simultaneously assess the total richness and the 
compactness of the cluster. 
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Lop - Lx{ROSAT) relation for REFLEX-NORAS clusters only 
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Table 7. The table presents the results of the best fit for the i and z band in the region 
of the minimum scatter, 0.2 < r < 1.0 Mpc h^Q^ for the Lop — Lx{ROSAT) relation in 
the subsample of REFLEX and NORAS X-ray selected clusters. We show the results for 
the orthogonal distance regression (ODRPACK) method. The a and /? parameters are 
given for several apertures with the 95% confidence errors. The orthogonal scatter and 
the scatters in log(ijf (0.1-2.4 keV)) and in log(Lop) to the best fit line are given by a, 
cflx and cfl^^ respectively. 

Finally by using the relations obtained from the z band, we demonstrated that, given 
the optical properties of a cluster, we can predict the X-ray luminosity and temperature 
with an accuracy of 60% and vice versa. By restricting the correlation analysis to the 
subsample of X-ray detected REFLEX-NORAS clusters, the minimum scatter decreases 
to less than 40% for the Lop — Lx relation. Since the observational uncertainties in the 
optical and in the X-ray luminosity are about 20%, the observed scatters in both relations 
should be intrinsic. 

The resulting logarithmic slope for the Lop — Lx relation with the minimum scatter 
is 0.38 ± 0.02, while the value for the Lop — Tx relation is 1.12 ± 0.08. Both results are 
not consistent with the assumption of hydrostatic equilibrium and a constant M/L. If we 
assume that M/L depends on the luminosity with the power law M/L oc L" '^ (Girardi 
et al. 2002), our results are in very good agreement with the expected values under the 
assumption of hydrostatic equilibrium. 

The analysis carried out in this paper on the correlation between X-ray and optical 
apparence of galaxy clusters is completely empirical. In principle, the best way to proceed 
in this kind of study is to measure the optical luminosity within the physical size of the 
cluster, like the virial radius. Without optical spectroscopic data or accurate temperature 
measurements, the chister virial radius can be calculated by assuming a theoretical model 
relating the optical luminosity to the cluster mass. At this stage of the work, we preferred 
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to tackle the cluster X-ray-optical connection with the empirical method explained in the 
paper, in order to have model-independent results. On the other hand, not taking into 
account the different cluster sizes could have affected both the slope and the scatter 
of the given relations (ea. If 41 IT^ and [T?)|l . Therefore, for a better understanding of the 
important connection between the X-ray and optical apparence of galaxy clusters, the 
optical luminosity has to be calculated within the physical size of the cluster. This work 
is in progress and will be published in the second paper of this series about the RASS- 
SDSS galaxy cluster sample. The next step will be the study of the fundamental plane of 
galaxy clusters. Through this kind of analysis we will find out if the observed scatter in 
the correlations between the optical and X-ray properties depends on another parameter 
related to the cluster compactness. Moreover, because of the link between the galaxy 
cluster fundamental plane and the M/L parameter, we will connect directly the slope of 
two relations to the behavior of M/L. 
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Wavelength (K.) 



Fig. 1. Response function of the SDSS photometric system. Dashed curves indicate the 
response function including atmospheric transmission at 1.2 airmass at the altitude of 
Apache Point Observatory. 
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Fig. 2. The plot shows a cluster region with the local background. The dots represent the 
galaxies in the sample. The biggest dots correspond to the brightest galaxies in apparent r 
magnitude. The local background number counts have been calculated inside the annulus 
with inner radius equal to r2oo + 0.2 deg and a width of 0.5 degree. The regions with voids 
due to lack of data or close clusters have to be discarded in the background estimation. 
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r 

Fig. 3. Global background number counts as a function of magnitude in the r band. The 
error bars include the contribution of large scale structure. The line shows the counts- 
magnitude relation expected for a homogeneous galaxy distribution in a universe with Eu- 
clidean geometry: N{r) = ArlO^-^^'^-'^^K The value of = 5.99±0.52{0.5mag)-^deg-'^ 
is the results of the fit in Yasuda et al. 2001 for 12 < r < 17. 
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Fig. 4. Example of the individual fitted luminosity function of a cluster in 4 Sloan 
photometric bands. The value of the fitted parameters are indicated in each panel. 
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Fig. 5. Comparison of the a parameter (slope) of tlie individual Schechter luminosity 
functions, aj, is calculated by using the local background correction, while ac is the 
result of the global background correction. The different background subtractions give 
consistent results. The error bars in the plot are at 68% confidence level. 
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Fig. 6. Comparison of tlie M* parameter (knee) of the individual Schechter luminosity 
functions. M*l is calculated by using the local background correction, while M*g is 
the result of the global background correction. As in the case of the slope, the different 
background subtractions give consistent results. The error bars in the plot are at 68% 
confidence level. 
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Fig. 7. Fit parameters of the individual Schechter luminosity functions in the r band. 
The fitting procedure is performed by using the local background magnitude number 
counts. The parameters show a slight correlation. The error bars in the plot are at 68% 
confidence level. 
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Fig. 8. Distribution of the reduced of the individual Schechter luminosity functions. 
Since the MLM allows us to perform a fit without binning the data but does not give 
information about the goodness of the fit, we have performed the fitting procedure with 
MLM and checked the goodness of our fitted luminosity functions with a test. All the 
cluster are well represented by the individual fitted luminosity function. 
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Fig. 9. The plot shows the correction for incompleteness to the total optical luminos- 
ity. The completeness magnitude limit is 21 mag in each Sloan photometric band. The 
correction is of the order of 5% for 50% of the cluster in the sample, while it is less than 
10% for 85% of the sample. The trend in the plot is due to a selection effect, since the 
most distant cluster are also the most luminous ones. 
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Fig. 10. Comparison of tlie cluster optical luminosities calculated from the cluster mag- 
nitude number counts with different background corrections. Ll is calculated with the 
local background subtraction, while Lq with the global one. The different corrections 
do not affect the cluster luminosity estimation. The error bars in the plot are at 68% 
confidence level. 
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Fig. 11. Comparison of the cluster luminosities determined with count-based and fit- 
based method. The error bars in the plot are at 68% confidence level. 



40 



P. Popesso et al.: ROSAT-SDSS Galaxy Clusters Survey. 



T3 1.4 
CD 

CO 
I 

W 1.2 

o 
o 

1 



m 

CO 0.8 



1 1 1 — 1 — I I I I I 

_r band 



I 



0.6 - 



"T 1 — I — I — I I I I I 1 1 — 1 — I — I I I I I 1 1 r 



lO'" 



10" 



10« 



10'3 



L (counts-based) (L^) 



Fig. 12. Comparison of the cluster luminosities count-based and fit-based. The plot 
shows the ratio of the two luminosities (count-based and fit-based) versus the count-based 
one. For 70% of the sample the ratio L fn-based/ L count-based is less the 1.1, indicating 
that fit-based luminosity is systematically less bright the count-based one. The reason is 
due to the Bright Cluster Galaxies (BCG): they are included in the computation of the 
count-based luminosity but not in the estimation of the fit-based one. 
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Fig. 13. The four panels show the behavior of the fit parameters in the King profile fit 
procedure. The first three panels show that there is no correlation between the parame- 
ters, while the bottom right panel presents the histogram of the exponent (3 of the King 
profile. 
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Fig. 14. The physical total size (the cross) of each cluster is estimated as the radius 
where the galaxy number density within the cluster becomes 3 times the error in the 
statistical background galaxy number density (the dashed line). 
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Fig. 15. Orthogonal scatter to the best fit of the Lopt — Lx relation obtained with the 
orthogonal method, as a function of the cluster aperture. The scatter shows a minimum 
in the region 0.2 < r < O.SMpc hyg^. In the plot different symbols are related to the 
different photometric bands in which Lgp is calculated: squares for the g band, triangles 
for the r band, hexagons for i and filled circles for the z band. The i and z bands show 
clearly the smallest scatter at any aperture. 
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Fig. 16. Orthogonal scatter to the best fit of the Lopt — Tx relation obtained with the 
orthogonal method, as a function of the cluster aperture. The scatter shows a minimum 
in the region 0.2 < r < 0.8 Mpc h^Q^. In the plot different symbols are related to the 
different photometric bands in which Lgp is calculated: squares for the g band, triangles 
for the r band, hexagons for i and filled circles for the z band. The i and z bands show 
clearly the smallest scatter at any aperture. 
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Fig. 17. Orthogonal scatter to the best fit of the Lopt — Lx relation obtained with the 
orthogonal method, as a function of the cluster aperture for the subsample of REFLEX- 
NORAS clusters. The scatter shows again a minimum in the region 0.2 < r < O.SMpc 
h^Q . In the plot different symbols are related to the different photometric bands in which 
Lop is calculated: squares for the g band, triangles for the r band, hexagons for i and 
filled circles for the z band. The i and z bands show clearly the smallest scatter at any 
aperture. 




Fig. 18. Correlation between optical luminosities and X-ray luminosities. The fit is 
performed with a linear regression in the log{Lopt) — log{Lx) space for each of the 4 
optical bands. The solid and the dashed lines are the results of the orthogonal and bisector 
method respectively over the all RASS-SDSS galaxy cluster sample. The dot-dashed line 
is the best fit result of the orthogonal method applied to the subsample of strictly X-ray 
selected REFLEX-NORAS clusters. The error bars are at the 68% confidence level in 
both the variables. 
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Fig. 19. Correlation between optical luminosities and ICM temperature. The fit is per- 
formed with a linear regression in the log{Lopt — log{Tx)) space for each of the 4 optical 
bands. The solid and the dashed lines are the results of the orthogonal and bisector 
method respectively. The error bars are at the 68% confidence level in both the vari- 
ables. 



